Quantum- "classical" correspondence in a nonadiabatic-transition system 
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A nonadiabatic-transition system which exhibits "quantum chaotic" behavior [Phys. Rev. E 63, 
066221 (2001)] is investigated from quasi-classical aspects. Since such a system does not have a 
naive classical limit, we take the mapping approach by Stock and Thoss [Phys. Rev. Lett. 78, 578 
(1997)] to represent the quasi-classical dynamics of the system. We numerically show that there is 
a sound correspondence between the quantum chaos and classical chaos for the system. 
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Nonadiabatic transition (NT) is a very fundamental 
concept in physics and chemistry 0,0- In atomic, molec- 
ular, and chemical physics literature, NTs occur as a 
breakdown of the Born-Oppenhcimer (BO) approxima- 
tion, which is essentially an adiabatic approximation to 
solve quantum systems with many degrees of freedom. It 
is still a tough problem to analyze the properties of NTs 
especially for multidimensional systems. One reason pre- 
venting us from deeper understanding of NT is the lack 
of a naive quantum-classical correspondence for NT like 
tunneling phenomena 0, Q . One way to get the classical 
picture for a NT system is to go back to the original 
system before the BO approximation: Fuchigami and 
Someda investigated dynamical properties of H^ from 
classical points of view by treating an electron and nu- 
clei as dynamical variables p|. Though there is a full 
uantum study for such a small system to compare with 
this "purist" way cannot be easily applied to much 
more "complex" systems. 

The mapping method recently advocated by Stock and 
Thoss can circumvent this deficiency. (This is reminis- 
cent of the Meyer-Miller method ||.) Their method is as 
follows: After the BO approximation, the discrete elec- 
tronic degrees of freedom are mapped onto the Schwinger 
bosons |3 • Since all the degrees of freedom become just 
bosons, the total system is rather easily treated semiclas- 
sically or quasi-classically. Using this method semiclas- 
sically, one can obtain, e.g., absorption spectra even for 
a pyrazin molecule with 24 degrees of freedom 10]. One 
can use it quasi-classically by solving the equations of 
motion derived from a mapping Hamiltonian. This is a 
very easy way to simulate NT systems because the addi- 
tional number of degrees of freedom for electronic parts 
is rather small. Using the periodic orbit theory [IlJ or 
the adiabatic switching method 01 ■> ° ne can obtain even 
quantum eigenenergies and eigenstates, in principle [l3j . 

On the other hand, multidimensional NT systems like 
Jahn- Teller molecules |l l| are known to show "quantum 
chaotic" behavior Fujisaki and Takatsuka investi- 

gated this problem deeply employing the two-mode-two- 
statc (TMTS) system which is considered as a system 



with two vibrational modes and two electronic states |15| . 
They calculated the statistical properties of the eigenen- 
ergies and eigenfunction for the TMTS system, and found 
that the system becomes strongly "quantum chaotic" un- 
der a certain condition. In addition, they showed that 
the chaos is not just a reflection of the lower adiabatic 
system nor that of the diabatic systems. (On the other 
hand, the chaos of all previous studies is just a reflec- 
tion of the lower adiabatic systems |14j.) This means 
that conventional classical descriptions do not help to ex- 
plain the quantum chaotic behavior. Hence this system 
deserves to be further studied from the "mapping" (ex- 
tended classical) points of view. Though there are some 
studies which investigated chaotic properties of this kind 
of mixed quantum-classical systems 16], our focus here 
is a quantum- "classical" correspondence (if any) for the 
TMTS system. 

The TMTS system first introduced by Heller [l7| 
is described by the following Hamiltonian: 



Htmts — 



Tkin + Va 
J 



J 



kin 



V B 



(1) 



where Tki n is the kinetic energy, Vi (i = A, B) is the 
potential energy for state i defined by 



V A 
V b 



(Pl+P 2 y), 

/ 2 2 , 2 2\ i 

(u x x +uj y y ) + 6.4, 
(u& 2 +uW) + e B , 



(2) 
(3) 
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with 
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£ = (x + 2a sin 0) cos 26> + (y - 2a cos 9) sin26>, (5) 
rj = -(x + 2a sin 0) sin 20 + (y - 2a cos 0) cos 20. (6) 

Note that we have just used a harmonic potential for 
each state. For the geometrical meaning of the param- 
eters, see Fig. ^ Here the Duschinsky angle and 
the nonadiabatic coupling constant J are two important 
parameters for the system. If these are appropriately 
chosen, the system becomes strongly "quantum chaotic" , 
i.e., the nearest-neighbor energy-level spacing distribu- 
tion becomes the Wigner type and the amplitude distri- 
bution of the eigenstates becomes Gaussian . 
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FIG. 1: A shematic representation of the TMTS system. The 
distance between the minima of the potential is 2a, and the 
angle between the relevant crossing seam (dotted line) and 
the primary axis of each potential (dashed line) is 9. Inset: 
The perspective of the TMTS system. The potential minima 
are different with Ae = es — e.4 = 0.173. 



with 

Ad(t) 2 = Ax A (t) 2 + Ap A (t) 2 + Ax B (t) 2 + Ap B {t) 2 
+ Ax(t) 2 + Ap x (t) 2 + Ay{t) 2 + Ap y (t) 2 (16) 

where x A = x A /y/2 + 2 7- Pa = Pa/V 2 + 2 7> x B = 
Pb = p A /^/2^, x = uj x x/V2E, y = 
uj y y/V2E, Px = Px/V2E, P y = p y /V2E, and E is 
the total energy for the system. Here tilde variables are 
introduced for normalization and A means a distance 
between a trajectory and its auxiliary one. In the fol- 
lowing we take the typical value of 7, i.e., 7 = 1/2 0. 
For the numerical calculation of Lyapunov exponents, the 
method by Benettin et al. |2Gj | is employed, i.e., we cal- 
culate A max (? 1 ) for a finite T, then shorten the distance 
Ad(T) to Ad(0), and run the trajectory again and so on. 
In this study, we took T — 24 because of the numerical 
divergence of the Ad(T) for larger values of T. 



The mapping Hamiltonian 0, for this system is 



^map 

with 



Tkin + N a Va + N B V B + J{x A x B + p A p B ) (7) 



Ni 



■Pi 



-1 (i = A,B) 



(8) 



Note that N A + Nb = 1. From this relation, one might 
consider Ni (i = A, B) as a probability but it is not 
the case. This is because the numerical range for Ni is 
—7 < Ni < I + 7. The equations of motion for this 
Hamiltonian are derived as follows: 



dt XA 



dt XB 



dt 
d 

~dt V 



V A p A 
VbPb 



Jpb, -r x p A = -V A x A - Jx B , (9) 
dt 



JPA, 



-rPB = -VbXb - Jx A (10) 
dt 



Py, dt Pv 



dVA N A -^N B , 



dx 
dV A 



N A 



dx 
dV B 



(11) 
(12) 



dy dy 

As a reference system, we take the lower adiabatic sys- 
tem as in [15| . The lower adiabatic systems is defined by 
the following Hamiltonian: 

(13) 



H. 



ad 



kin 



v. 



ad ' 



with 



Va + Vb 



'V a -Vb 



.. 1 



J 2 



(14) 



Since the TMTS system can be quantum chaotic as 
mentioned above, we investigate the chaotic properties 
for the mapping system and the lower adiabatic system. 
As an indicater of chaos, it is very natural to take Lya- 
punov exponents [l9l |. We focus on a finite-time maxi- 
mum Lyapunov exponent calculated as 



ccn^ioj 



Ad{T) 
Ad(0) 



(15) 
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FIG. 2: J dependence of the distribution of the Lyapunov 
exponents for the mapping system. The nonadiabatic cou- 
pling is (a) J = 0.3, (b) J = 1.5, and (c) J = 7.5. The finite 
time is T = 24, and the iterative number of the time average 
is 10. The Duschinsky angle is 9 = 7r/3. 
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FIG. 3: J dependence of the distribution of the Lyapunov 
exponents for the lower adiabatic system. The nonadiabatic 
coupling is (a) J = 0.3, (b) J = 1.5, and (c) J = 7.5. The 
finite time is T = 24, and the iterative number of the time 
average is 10. The Duschinsky angle is 9 — ir/3. 



FIG. 4: 9 dependence of the phase space averaged Lyapunov 
exponent for the mapping system with J — 1.5. The Duschin- 
sky angle is (a) 9 = 0.0, (b) 9 = tt/6, and (c) 9 = tt/3. The 
finite time is T = 24, and the iterative number of the time 
average is 10. 



Following the previous studies |l5j, we concentrate on 
a rather high energy region around E = Eq = 28. (This 
energy is much higher than that around the crossing seam 
region.) Varying the Duschinsky angle 8 and the nonadi- 
abatic coupling constant J, we calculate the distribution 
of the finite-time maximum Lyapunov exponent A max (T). 
For simplicity, hereafter, we call A max (T) just as a Lya- 
punov exponent. Since the phase space has a structure 
especially for the lower adiabatic system, globally aver- 
aged Lyapunov exponents are not so useful. Instead, 
we investigate the properties of the distribution of the 
Lyapunov exponents, which reflects phase space struc- 
ture of the system. (Remember that the Berry-Robnik 
distribution, which reflects the phase space volume of 
chaotic seas, is useful for mixed systems |2l|.t To this 
end, we prepare an initial ensemble of particles which 
are sampled from a part of the equi-energy potential sur- 
face Va{x,i/) — Eq with constraints p x = p y = and 
y < xtan# + a/ cos 8. (The latter constraint means that 



we only take points below the crossing seam. See Fig.[3) 
We take 40 sample points from this curve and calculate 
the histogram for the Lyapunov exponents. We believe 
that these sample points represent a typical situation of 
the TMTS system because, at least, the characteristic of 
the lower adiabatic system can be understood from these 
sampling points |15|. 

First we fix 8 — tt/3, and investigate the J dependence 
of the chaotic properties. Let us summarize the corre- 
sponding quantum results: the nearest-neighbor (energy- 
level) spacing distribution is Wigner for J = 1.5 (3rd 
row of Fig. 3 in whereas it is rather Poisson for 

J = 7.5 (3rd row of Fig. 4 in Ha h and rather mixed for 
J = 0.3 (3rd row of Fig. 2 in pjl. Figure H shows the 
distribution of the Lyapunov exponents for the TMTS 
system. The distribution for J — 1.5 has a sharp peak 
around 1, whereas that for J — 7.5 has a rather broad 
peak around 0.4, and that for J = 0.3 is also a little 
bit broad. This corresponds to the quantum results, at 
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least, qualitatively. On the other hand, Fig. [3] shows the 
distribution of the Lyapunov exponents for the lower adi- 
abatic system. As one can see, the values themselves are 
much smaller than those for the TMTS system and it is 
difficult to distinguish three distributions. This also cor- 
responds to the quantum mechanical calculation for the 
lower adiabatic system (Fig. 9 in 0]). From this com- 
parison, it is reasonable to conclude that there is a sound 
quantum- "classical" correspondence between the TMTS 
system and the mapped system in view of their "chaotic" 
properties. 

Next we investigate the 9 dependence of the chaotic 
properties while fixing J = 1.5. Although there is a 
strong peak around A max (T) — 1.3 as shown in Fig. 0] 
(a), we can see that the system with 9 = is not globally 
chaotic. This is because x does not effectively coupled 
to y when 9 = 0, and the motion along x axis is regular. 
(On the other hand, even with 9 = 0, the motion along y 
axis can be chaotic as shown by the peak of the distribu- 
tion around 1.3.) In such a case, we do not expect that 
the Wigner type distribution arises in the correspond- 
ing quantum system, and this is the case for the TMTS 
system (1st row of Fig. 3 in 0]). On the other hand, 
for intermediate Duschinsky angles (9 = 7r/6,7r/3), the 
Lyapunov exponent distributions show that the system 



is globally chaotic [Fig.0](b),(c)], and the corresponding 
quantum system can have the Wiger type distribution, 
which is also confirmed numerically (2nd and 3rd rows 
of Fig. 3 in ^f|). Of course, we have to admit that this 
correspondence is loosely stated, and there remains a dif- 
ficult question exactly when the quantum chaos begins in 
the parameter space. Such an issue must be addressed 
utilizing semiclassical methods 0, El or the adiabatic 
switching method ^ is also interesting to investi- 

gate phase space structure of this mapping system and 
to relate it to the nearest-neighbor spacing distribution 
via the Berry-Robnik distribution |2lJ . We also hope that 
this study will cast a light on the relation between the 
statistical reaction theory for NT systems and Lyapunov 
spectra for them. 

In this paper, employing the mapping approach 
by Stock and Thoss, we investigated a nonadiabatic- 
transition system which exhibits "quantum chaotic" be- 
havior from quasi-classical aspects. By comparing the 
statistical properties of the quantum system with the 
Lyapunov exponent distributions of the mapping system, 
we found that there is a sound quantum- "classical" cor- 
respondence in the system. 

The author thanks T. Takami for comments. 



[1] H. Nakamura, Nonadiabatic Transition: Concepts, Ba- 
sic Theories and Applications (World Scientific, Singa- 
pore, 2002); C. Zhu, Y. Teranishi, and H. Nakamura, 
Adv. Chem. Phys. 117, 127 (2001). 

[2] K. Takatsuka, Y. Arasaki, K. Wang, and V. McKoy, Fara- 
day Discuss. 115, 1 (2000); Y. Arasaki, K. Takatsuka, 
K. Wang, and V. McKoy, Phys. Rev. Lett. 90, 248303 
(2003). 

[3] H. Ushiyama and K. Takatsuka, Phys. Rev. E 53, 

115 (1996); K. Giese, H. Ushiyama, and O. Kiihn, 

Chem. Phys. Lett. 371, 681 (2003). 
[4] A. Shudo and K.S. Ikeda, Physica D 115, 234 (1998); 

T. Onishi, A. Shudo, K.S. Ikeda, and K. Takahashi 

Phys. Rev. E 64, 025201(R) (2001). 
[5] S. Fuchigami and K. Someda, J. Phys. Soc. Jpn. 72, 

1891 (2003). See also M. P. Strand and W. P. Reinhardt, 

J. Chem. Phys. 70, 3812 (1979). 
[6] I. Kawata, H. Kono, and Y. Fujimura, 

J. Chem. Phys. 110, 11152 (1999). 
[7] G. Stock and M. Thoss, Phys. Rev. Lett. 78, 578 (1997); 

M. Thoss and G. Stock, Phys. Rev. A 59, 64 (1999). 
[8] H.-D. Meyer and W. H. Miller, J. Chem. Phys. 70, 

3214 (1979); ibid. 71, 2156 (1979). See also X. Sun and 

W. H. Miller, ibid. 106, 6346 (1997). 
[9] J.J. Sakurai, Modern Quantum Mechanics (Addison- 

Wesley, Reading, MA, 1994). 
[10] M. Thoss, W.H. Miller, and G. Stock, 

J. Chem. Phys. 112, 10282 (2000). 
[11] M. C. Gutzwiller, Chaos in Classical and Quantum Me- 
chanics (Springer- Verlag, New York, 1990). 
[12] W. P. Reinhardt, Adv. Chem. Phys. 73, 925 (1990); 

R.T. Skodje and JR. Cary, Comp. Phys. Rep. 8, 223 

(1988). 



[13] J. Zwanziger, ER. Grant, and G.S. Ezra, 
J. Chem. Phys. 85, 2089 (1986); M. Ple- 
tyukhov, Ch. Amann, M. Mehta, and M. Brack, 
Phys. Rev. Lett. 89, 116601 (2002); B. Balzer, S. Dilthey, 
S. Hahn, M. Thoss, and G. Stock, J. Chem. Phys. 119, 
4204 (2003). 

[14] H. Koppel, W. Domcke, and L.S. Cederbaum, 
Adv. Chem. Phys. 57, 59 (1984); D.M. Leitner, 
H. Koppel, and L.S. Cederbaum, J. Chem. Phys. 104, 
434 (1996); H. Yamasaki, Y. Natsume, A. Terai, and 
K. Nakamura, Phys. Rev. E 68, 046201 (2003). 

[15] H. Fujisaki and K. Takatsuka, Phys. Rev. E 63, 
066221 (2001). See also H. Fujisaki and K. Takatsuka, 
J. Chem. Phys. 114, 3497 (2001); H. Higuchi and 
K. Takatsuka, Phys. Rev. E 66, 035203(R) (2002). 

[16] R. Blumel and B. Esser, Phys. Rev. Lett. 72, 3658 
(1994); A. Bulgac and D. Kusnezov, Chaos, Solitions, 
and Fractals, 5, 1051 (1995); J. Gong and P. Brumer, 
Phys. Rev. Lett. 88, 203001 (2002). 

[17] E.J. Heller, J. Chem. Phys. 92, 1718 (1990). 

[18] See J. Tang, M.T. Lee, and S.H. Lin, J. Chem. Phys. 119, 
7188 (2003), and references therein. 

[19] E. Ott, Chaos in Dynamical Systems, 2nd edition (Cam- 
bridge Univ. Press, Cambridge, 2002). 

[20] G. Benettin, L. Galgani, and J. Strelcyn, Phys. Rev. A 
14, 2338 (1976). For recent accurate methods to 
calculate Lyapunov exponents, see T. Okushima, 
Phys. Rev. Lett. 91, 254101 (2003). 

[21] M.V. Berry and M. Robnik, J. Phys. A 17, 2413 (1984); 
H. Makino, T. Harayama, and Y. Aizawa, Phys. Rev. E 
63, 056203 (2001). 



